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Abstract 

The review chapter starts by a pedagogical introduction to the general concept of the scattering 
theory: from the fundamental wave-function picture to the second-quantization language, with the 
aim to clear possible ambiguity in conventional textbooks. Recent progress in applying the method 
to current fluctuations and oscillating-parameter driven quantum pumping processes is presented 
with inclusion of contributions by Biittiker, Brouwer, Moskalets, Zhu, etc. In particular, the spin- 
orbit-coupling affected shot noise can be dealt with by taking into account the spin-dependent 
scattering processes. A large shot noise suppression with the Fano factor below 0.5 observed 
experimentally can be illustrated by effective repulsion between electrons with antiparallel spin 
induced by the Dresselhaus spin-orbit coupling effect. A Floquet scattering theory for quantum- 
mechanical pumping in mesoscopic conductors is developed by Moskalets et al., which gives a 
general picture of quantum pumping phenomenon, from adiabatic to non-adiabatic and from weak 
pumping to strong pumping. 
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I. PEDAGOGICAL INTRODUCTION TO THE GENERAL CONCEPT OF THE 
SCATTERING THEORY 



A. Landauer-Biittiker conductance 



The discussion is based on the scattering approach to electrical conductance. This ap- 
proach, as we will show, is conceptually simple and transparent. Nevertheless, the generality 
of the scattering approach and its conceptual clarity, make it the desired starting point of 
a discussion of noise in electrical conductors. By expanding the time-dependent scattering 
matrix into Fourier series, a description of the quantum pumping phenomenon can be given. 

We start with the wave function picture and consider an electron tunneling through a 
one-dimensional single barrier 1 , which can be realized in a semiconductor heterostructure 
with a layer of Al a .Gai„ a; As of width 2L imbedded in GaAs as shown in Fig. 1. In the 
effective mass approximation, the electron motion in each layer of the structure is described 
by the stationary solution of the envelop equation in the x-direction. 

'-h 2 d 1 d 



+ Veff (X) 



ip (x) = Eip (x) . 



(1) 



2 dx m* (x) dx 

Here, m* and V e ff are the effective mass and potential in different regions with E the energy 
of the transporting electron. The electron's wave functions are expressible as 



^ (x,t) = ip (x) e 



-iEt/h 



(Ae tkx + Be" ikx ) e ~ iEt/r \ x < —L, 
(Ce KX + De KX ) e - iEtlr \ -L<x <L 
(Ee ikx + Fe~ ikx ) e - lEt/h , x > L, 



(2) 



with k = \Jlm*E j h and k = ^2m* (Vq — E) j h. The coefficients A and B are associated 
respectively with incoming and outgoing waves on the left side relative to the barrier. Like- 
wise, the coefficients E and F are respectively outgoing and incoming waves on the right. 
The scattering matrix connects the incoming and outgoing fluxes as 



(3) 



Ideal (i.e., without scattering) conducting leads connect the scattering region to reservoirs 
on the left and right characterized by quasi- Fermi energies fii and /z 2 , respectively, corre- 
sponding to the electron densities there. These reservoirs or contacts randomize the phase of 
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the injected and absorbed electrons through inelastic processes such that there is no phase 
relation between particles. For such an ideal ID system, the current injected from the left 
and right may be written as an integral over the flux->2 
2e 



2tt 



dkv (k) h (k) T (E) 



dk'v (k') f 2 (k') T (E') 



(4) 



where v(k) is the velocity, T(E) is the transmission coefficient, which can be obtained from 
the wave function Eq. (2), and f\ and f 2 are the reservoir distribution functions characterized 
by \i\ and fi 2 , respectively. The integrations are only over positive k and k! relative to the 
direction of the injected charge as positive k is in +x-direction and positive k' is in — re- 
direction. If we now assume low temperatures, electrons are injected up to an energy \i\ 
from the left lead and injected up to \ii from the right one. Converting to integrals over 
energy, the current becomes 

II = I [JT dE (§) v (k) T (E) - f» dE (%) v (k 1 ) T (E)] 
= ^.i:;dET(E). 

It can be seen that the first term of Eq. (5) is the flux generated by electrons injected from 
the left reservoir and the second term is that from the right reservoir. The integration is 
done in two independent ensembles. 

The Landauer-Btittiker conductance introduced above can be reproduced in the second- 
quantization language. Without loss of generality, we assume the input amplitudes from the 
two reservoirs in Eq. (2) to be unity. For simplicity, as a single channel is considered, the 
scattering matrix can be described in the relation 
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or elaborately 
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where operators <2l/r annihilate electrons incident upon the sample from the left /right reser- 
voir and operators bL/R describe electrons in the outgoing states, t and t' are the transmis- 
sion amplitudes of electrons incident rightward and rightward, respectively, and r and r' are 
the corresponding reflection amplitudes, as defined conventionally. Hence, the Landauer- 
Biittiker formula of the current can be expressed as 

e 



2ttH 



dE{ 



i{ (E) a L (E) - b[ (E) b L (E) 



(8) 



where (• • • ) calculates the quantum statistical average of the product of an electron creation 
operator and annihilation operator of a Fermi gas. As a conventional electron reservoir is 
considered, we have 



4 (E) a L (E')) = f L (E)5(E-E'), 

a R (E)a R (E>)) = f R (E)5(E-E>), (9) 
4 (E)a R (E')) = (a R (E)a L (E')) = 0. 



It should be noted from the third formula of Eq. (9) that electrons incident from the 
left and right reservoirs are completely incoherent and the phases of the two reservoirs are 
randomized and completely unrelated. Substituting Eq. (9) to Eq. (8), the formula of the 
current expressed in Eq. (5) can be reproduced. 

It is interesting to see that if we could build up a system with the conductor connected 
to two correlated reservoirs, in which the quantum statistical average of the cross product 
reads 

'4 (E) a R (£')) = (4 (E) a L (£')) = fs (E) 5 (E - E>) , 

1, fi < min (/ii,/i 2 ) , ( 10 ) 
0, others, 

the system can generate a current demonstrating the interference between the electron states 
incident from the left reservoir and the right one. 
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(11) 



In the next subsection, we would use a toy system based on correlated reservoirs to further 
illustrate the scattering scheme. Significant difference between uncorrected and correlated 
reservoirs is illuminated. It is also elaborated that in scattering problems the transmission 
and reflection process demonstrates the quantum state spanned throughout the space. 
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B. Further illustration of the scattering scheme in a toy system with correlated 
reservoirs 



We consider incidence only from the left in Eq. (2) and the input flux is assumed to be 
unity. The wave function in different scattering regions can be expressed as 



i) (x) = < 



e ikx _|_ re ikx ^ x < 

Ce KX + De~ KX , -L<x<L, (12) 
te lkx , x > L. 



2 1 2 

An incident particle is transmitted with probability \t\ and reflected with probability |r| . 
It is determined by the wave function that the particle momentum has value of hk with 
probability (l + |t| 2 )/2 and —hk with probability |r| 2 /2. The mean value of the momentum 
\t\ hk characterizes the density of the probability current J, which can also be obtained 
from the continuity equation dp/dt + V • J = with p the probability density. 

An interesting interference pattern can be observed if we consider the incidence from the 
left and the right is correlated. We consider a toy system of a one-dimensional electron 
gas subject to two oscillating gate volatages (see Fig. 2). The inspiration comes from the 
quantum pumping phenomenon. The single-particle Hamiltonian reads 

H= -£w +u (13) 

with U (x, t) = Q{x + 2L) (-L - x) U x (t) + (x - L) (2L - x) U 2 (f) . The two barriers 
Ui and U 2 are adiabatically modulated at the frequency u> with a phase difference 0. 

Ui(t) = U 10 + UiuWiu;t, ^ 
U 2 (t) = U 20 + U 2w sin (ut -(f)). 

We assume u is extremely small so that a static treatment is valid. Therefore, it is tolerable 
to consider only the zero order of the Fourier component of the time-dependent scattering 
matrix without taking into account the photon- absorption/emmission processes, which is 
done in adiabatic quantum pumping theory. 

As shown in Fig. 2, the electrons are incident from the left and right reservoirs with 
identical amplitudes at zero bias, which is set to be unity without impairing generality. As 
assumed, incidence from the two correlated reservoirs characterize a coherent single-particle 
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state. The single-particle wave function at a certain time has the following form. 



< 



[e ikx + (r + t'e id ) e' ikx ] e ~ iEt / h , x < -2L, 

[A 2 e K2X + B 2 e~ K2X ] e~ iEt / r \ -2L <x< —L , 

[A 3 e ikx + B 3 e~ ikx ] e~ iEt / h , -L < x < L , (15) 



[A 4 e K4X + B 4 e~ K4X ] e " iEt ' r \ L < x < 2L , 
[(t + r'e id ) e ikx + e~ ikx ] e - iEt ' r \ x > 2L. 



k = \ / 2m*E J h and K2/4 



2m* (U\/2 — E^j h. t and r quantify the transmission and 



reflection amplitudes of the electrons incident from the left reservoir while t' and r' quantify 
those incident from the right with i! = t and r' = —tr*/t* . We introduce a geometric 
phase 6 to describe the unavoidable phase difference between the electrons injected from 
the two reservoirs. In the toy approach, correlated reservoirs are assumed, which justifies a 
particular value of 9. For simplicity and without violation of the physics we take 9 = 0. It 
is noted that in conventional real reservoirs mixed-ensemble integral should be applied, i.e., 
the probability flow of the electron incident from the left reservoir absolutely cancels out 
that of the incidence backward at zero bias as a result of the randomized phase distribution, 
which absolutely differs from our toy consideration. 
From the continuity equation 



we could derive the probability current flow as functions of the transmission and reflection 
amplitudes as 



We can also see from the wave function Eq. (15) that the particle momentum has value of 
hk with probability [1 — Re (tr*)]/2 and — hk with probability [1 + Re (tr*)]/2. The mean 
value of the momentum —Re (tr*) hk characterizes the density of the probability current of 
Eq. (17). The net current density can be described by the period- average of the probability 
current density multiplied by the carrier charge and density. The accumulated contribution 
by electrons within the ±hu sidebands is taken into account by an integral. Without dynamic 
modulation, no current occurs as the ±hu energy channel is closed even when the probability 
current density is nonzero for asymmetric barrier configuration. The current density as a 




(16) 




(17) 
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function of the Fermi level thus becomes 



h {E F ) 



Ah f Ep+hw ekp e N (E) f (E) 
™* J Ep—hw (27r/w) 



L 



Re [r* (t) t' (t)} dtdE, 



(18) 



where the density of states of a one- dimensional electron gas is 



N(E) 



V jm e 

'nhy^E' 



(19) 



Here, k is the wave vector of the electron, e is the electron charge. p e is the carrier density of 
the two-dimensional electron gas in which the quantum wire is confined. f(E) is the Fermi- 
Dirac distribution function of the leads. V quantifies the volume of the one- dimensional 
wire and m e is the mass of a free electron. 

Here, in the toy configuration, quantum interference is remarkably demonstrated. The 
single-electron state Eq. (15) interferes with itself and carries the probability flow and hence 
the net current through oscillating cycles. 

We numerically calculated the current in a one-dimensional electron gas (2DEG) based 
on the GaAs/AlGaAs heterostructures with the average carrier density p e ~ 10 11 cm -2 and 
the effective mass of the electron m* ~ 0.067me ? . The width of the two gate potential 
barriers L = 20 A equally separated by a 2L = 40 A width well. The amplitudes of the 
modulations Ui u = = 1 meV. All of the above setups are not essential as we are dealing 
with an assumed toy structure. 

In Fig. 3, it is shown that the time-integrated current demonstrates a sinusoidal pattern 
as a function of the phase difference between two oscillating parameters, which is a result 
of quantum interference. The probability density flow formulated in Eq. (17) is a result of 
phase difference between transmission forward and backward. Quantum phase interference 
gives rise to nonzero probability density flow for asymmetric barrier configurations. The 
absolute strength of the probability flow increases as the height difference between the two 
potential barriers increases determined by the phase difference 0. Fig. 4 presented the 
time variation of the net current within an oscillating cycle. Considering the time-averaged 
effect, the integrated asymmetry is different. When the phase difference between the two 
modulations approaches tt, time-reversal symmetry destroys the time-integrated current to 
zero although the probability density flow maximizes at a certain time. i.e. the probability 
flow in half a pumping cycle completely offsets that of the other. Therefore a sinusoidal 
dependence on <p occurs in the time-integrated current. 
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In real reservoirs, the single-particle state between reservoirs has a definite momentum 
direction determined by its source. When coherent reservoirs can be realized in any form, 
however, quantum states within the mesoscopic conductor can be expressed as Eq. (15) and 
double-slit interference pattern is observable in an electron device. 

C. Time-dependent scattering- matrix theory 

The scattering-matrix equation b a = S a pap with a and (3 indexes of lead, channel, and 
spin introduced in Sec. LA. for the Landauer-Biittiker conductance characterizes the trans- 
port properties through a conductor at a certain bias. 

In a more general situation with dynamic processes, e.g. in quantum pumping, a time- 
dependent scattering matrix can be introduced as follows. 

/oo 
S a p(t,t')ap(t)dt', t>t, (20) 
-oo 

with a and f3 general indexes denoting the lead, channel, and spin. An incident state dp at 
time t' is scattered into the outgoing state b a at time t with the amplitude S a p (t,t'). 

The time-dependent scattering-matrix picture described by Eq. (20) is exactly equiva- 
lent to the time-dependent Schrodinger equation with the elements of the scatting matrix 
amplitudes of the wave function. In usual cases, the time-dependent Schrodinger equation 
cannot be solved exactly, similarly to the time-dependent scattering matrix. In the static or 
adiabatic cases, it is advantageous to use an analog of the Wigner transform for the matrix 

S a j3 (t, t')-, 

/oo 
e iE ^S aP (t,t')dt'. (21) 
-oo 

An on-time scattering process S a p (t, t') with t = t' is sufficient to describe the bias-driven 
conductance. Namely, from Eq. (21), we can use S a p (E) with E labeling the energy channel 
to fully capture the transport physics. 

When the scattering time t — t' is small (i.e., the dynamic characteristic frequency is 
much smaller than the inverse Wigner time delay), the dynamics can be approximated into 
the instant-scattering picture. Physically this means that the scattering matrix changes 
only a little while an electron is scattered by the mesoscopic sample under dynamic modu- 
lation, in which we use the term "adiabatic" . In adiabatic dynamics, we can use low-order 
Fourier components of S a @ (E, t) to characterize transport physics. Therefore, small t — t' is 
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transformed into variation of the particle energy by side-band broadening around the Fermi 
level. 

In Sec. Ill, we would illustrate the time-dependent scattering approach in adiabatic 
quantum pumping beyond the linear-response approximation. 

II. SPIN-ORBIT COUPLING AFFECTED SHOT NOISE 
A. Background 

Current fluctuations are present in almost all kinds of conductors and have been devel- 
oped into a very active and fascinating subfield of mesoscopic physics (for review see Ref. 



5|). At low temperatures, thermal fluctuations are extremely small, the current fluctuation 
properties are governed by the so-called shot noise, which is a consequence of the quantiza- 
tion of charge. Shot noise is useful to obtain information on a system which is not available 
through conductance measurements. In particular, shot noise experiments can determine 
the quantum correlation of electrons, the charge and statistics of the quasi-particles rele- 
vant for transport, and reveal information on the potential profile as well as internal energy 
scales of mesoscopic systems^ - — . Shot noise is generally more sensitive to the effects of 
electron-electron interactions than the average conductance. 

A convenient measure of shot noise is the Fano factor F, which is the ratio of the actual 
shot noise and the Poisson noise. The Poisson noise would be achieved in measurement if 
the transport is carried by single independent electrons. Four typical values of the Fano 
factor characterize the shot noise properties of different mesoscopic conductors. 

1. F = 1 characterizes Poissonian processes. Particles are completely independent during 
transport corresponding to channels through which transmission is exponentially small. In 
diffusive transport, they are the so-called approximately-closed channels. Typical conductors 
featuring F = 1 include tunneling junction, Schottky-barrier diode, and asymmetric double- 
barrier diode. 

2. F = characterizes ballistic transport. In ballistic transport, transmission approaches 
the maximum of unity. Free particles wave function extends throughout the space, i.e. 
particle beams exhibit full coherence. In diffusive transport, they are the so-called open 
channels. Typical conductors featuring F = include pure metal and free two-dimensional 
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electron gas. 

3. F — 1/3 characterizes diffusive transport. Open channels and closed channels are 
distributed randomly. As a result of ensemble average of channels, the strength shot noise 
falls between F — 1 and F — 0. Typical conductors featuring F = 1/3 include diffusive 
metals, graphene, and two-dimensional electron gas modulated by magnetic barriers. 

4. F = 1/2 characterizes ballistic transport constrained by Pauli principle. The Pauli 
principle forbids two electrons to be in the same channel simultaneously. As a result, the 
shot noise is suppressed. All kinds of conductors are subject to the Pauli principle. In 
the symmetric-double-barrier diode, where Pauli exclusion is the only correlation between 
particles, the shot noise features F = 1/2. 

With increased attention^— to semiconductor spintronics, materials such as GaSb, InAs, 
and InSb with considerably strong spin-orbit coupling (SOC) const ant^ 1 ^ are becoming 
widely used in mesoscopic conductors, observations^ 1 ^ beyond the formalism of earlier 
theory on current shot noise are reported. The scattering approach is developed to derive 
a general formula for the shot noise in the presence of the SOC effect and apply it to the 
double-barrier resonant diode (DBRD) systems. It is demonstrated that the microscopic 
origin of the super-suppression of the shot noise observed in experiments^ is the bunching 
interaction between electrons with opposite spins resulting from the Dresselhaus k 3 terms^ 1 ^ 
in the effective Hamiltonian of the bulk semiconductor of the barriers. 

B. Theoretical approach 

In this part, the effect of the Dresselhaus SOC to the shot noise properties in the DBRD 
structure connected to ferromagnetic or normal metal leads is considered. The theory can be 
generalized to other coherent mesoscopic conductors subject to Dresselhaus and/or Rashba 
SOC. 

The Dresselhaus SOC is caused by the bulk inversion asymmetry and exists broadly in 
III-V compound semiconductors with zinc-blende crystal structures^. Consider the trans- 
mission of electrons with identical wave vector k = (kn,A; z ) through a certain potential 
barrier grown along z || [001] direction, ky is the wave vector in the plane of the barrier 
and k z is the wave vector component normal to the barrier. The electron Hamiltonian of 
the barrier in the effective-mass approximation contains the spin-dependent k 3 Dresselhaus 
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term 

d_ 

dz 



H D = j(cr x k x - cryky)^, (22) 



where 7 is the material constant denoting the strength of the Dresselhaus SOC, a x and a y 
are the Pauli matrices. The Dresselhaus term can be diagonalized by the spinors 




X± = -z= > ( 23 ) 



which describe the spin-up ("+") and spin-down ("— ") electron eigenstates. 

Suppose the system is a layered mesoscopic conductor with its potential profile described 
by Vq(z) (see Fig. 5), the electron motion in each layer of the structure is described by the 
Hamiltonian 

H = -^-y + V(z)+H D , (24) 

where V(z) = Vq(z) — eF(z + b)Q(z + b)Q(a + c — z) with F the magnitude of the electric 
field, 0(2) the step function, and —h and a + c the longitudinal coordinates of surfaces in 
z direction. Our discussion is within the framework of single electron approximation and 
coherent tunneling^^, and only zero-frequency noise at zero temperature is considered. 
Under the assumption that k\\ is conserved during the tunneling, the wave functions for the 
electrons with definite longitudinal electron energy (E z ) can be obtained from Schrodinger 
equation based on the Hamiltonian given in Eq. (24), which can be diagonalized by spinors 
X±- So, the wave functions become 



cxp(/k| ■ ,>) ( Y^^^-exp{ik ljZ )xj + E r i\J^- e M-ihjz)xjJ , z<b, 

exp(zk|| • p) E $u( z )-> -b < z < a + c, 



exp(zk|| • p) E tjJ-g£:exp(ik Bj z)xj, z > a + c, 

j=± V 3 

(25) 

where $^ denotes the wave function in the conductor region, and tj and rj are the transmis- 
sion and reflection amplitudes, which can be calculated using the transfer-matrix method^. 
The spin "+" and spin "— " components of the electron wave functions transport separately 
without correlation. As standard scattering method is applied, we introduce creation and 
annihilation operators of electrons in the scattering states. Schematics of the scattering 
states are shown in Fig. 6. Operators o/ Lna (E) and a Lnu (E) create and annihilate electrons 
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with total energy E and spin polarization a in the transverse channel n in the left lead, 
which are incident upon the sample. In the same way, the creation b\ na {E) and annihilation 
bLna(E) operators describe electrons in the outgoing states. They obey anti-commutation 
relations. Therefore, we can write the scattering matrix of the sample as 

f r+ t' + \ 

r_ t'_ 

t+ r' + 

\ t- r'_ ) 
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M 2 



with *' = t,r' = -£r*, Tj (E) = \tj\ 2 , Rj(E) = |r,-| 2 , and j = ±. The matrices M 1 and 
M 2 are unitary transformations between spin "|4" states and spin "±" states, s is the 
scattering matrix connecting the incoming and outgoing spin "±" states of the nth channel. 
The current of the system can be derived as follows 



ap mn a' a" cr 



i(E-E')t/hrf 



ai m Amz^,AL-,E,E')d, n AE'), 



where 

k 



(27) 



(28) 



For a system at thermal equilibrium, the quantum statistical average of the product of an 
electron creation operator and annihilation operator of a Fermi gas with spin polarization is 

a Lt a ^"t) = hp, {^Ln^Lnij = J La, {^Ln^ a Lnij = {^Lni a Ln^ = \{fhp - J La) • (29) 

Without loss of generality, we set the unit vector directed along the spin orientation n s = 
(1, 0, 0). Making use of Eqs. (24)- (26), after some algebra, we can obtain the expression for 
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the zero-frequency noise power 
Sap = Safi(0) = £ EfdE 

662 ran a'cr" a 
CT2 

X { [Sa"a' 2 Sa'a^fs(E) + 5 CT »4<W 2 ' \ (fps(E) - fas(E))] 

- [S<r>*»fs(E) + <W»| (f pS (E) - f aS (E))] x [5 a/ , al J S2 (E) + (f P s 2 (E) - U(E))] } . 

(30) 

Eq. (30) gives a general formula to calculate the shot noise in mesoscopic conductors subject 
to the SOC effect. It can be used to predict the low-frequency noise properties of arbitrary 
mult i- channel, multi-probe phase-coherent conductors in the presence of Dresselhaus SOC. 
It can be naturally extended to the system with Rashba SOC and the system with both 
Dresselhaus SOC and Rashba SOC. In the limit of zero SOC and when the conductor is 
connected to normal metal leads, there is f aa (E) = f ap {E) and T n+ (E) = T n ^(E) = T(E), 
thus Eq. (30) reconverts to the formula provided by Biittiker 32 concerning scalar electron 
systems without the spin degree of freedom inducing observable difference. 

C. Numerical results in comparison with experiment 

To further demonstrate our theory, we provide numerical results based on real heterostruc- 
tures in comparison with experiment. 

The double-barrier structure considered here is constructed of layers of Ga x Ali_ x Sb with 
x = 0.15/0.3/0/0.3/0.15 (x = 0.3 for the barriers and x = for the well), which are 
known to be semiconductors with relatively strong Dresselhaus SOC^&Si. Our target setup 
is a symmetric double-barrier structure with the thickness of the well a = 30 A and the 
thickness of the two barriers b = c = 50 A. The height of the barrier = 230 meV 
and the depth of the well V w = 200 meV are given by the heterostructure properties^ (cf. 
Fig. 5). We assume that in the whole region the effective mass m* = 0.053771^^. The 
chemical potential of the two electrodes is set to be 12 meV. For comparison, we choose the 
Dresselhaus constant 7 = 0, 40, 80, 120 eV A 3 . 

We see the current fluctuation of the system caused by Dresselhaus SOC. Fig. 7 presents 
results of the electric current I and the shot noise S versus the external bias. It is demon- 
strated that the peaks of the current and of the shot noise are lowered as the Dresselhaus 
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constant 7 increases, and a concave down in the ascending side of the noise curve is obvious 
for non-zero 7. The concaveness gives rise to nadirs far below 0.5 in the Fano factor (see 
Fig. 8 (a)). To compare with experiment, S vs I curves normalized to unity are displayed in 
Fig. 8 (b). In the positive differential conductance region, the shot noise follows the value 
of uncorrelated electrons (2eJ) for small tunnel currents and is significantly suppressed for 
larger currents and eventually increases. The suppression is above one-half for 7 = 0, near 
one-half for 7 = 4.0 x 10~ 29 eV m 3 , and below one-half for 7 larger than 8.0 x 10 -29 eV m 3 . 
In the negative differential conductance (NDC) region, the Coulomb interaction or charging 
effect in the well enhances the shot noise and overweighs the effect of SOC 33 . Iannaccone et 
al. have focused on the NDC region and obtained enhanced shot noise^. 

The results shown in Figs. 7-8 can be understood from the following. The SOC interaction 
behaves like a pseudo magnetic field and induces split of different spin components of the 
resonant level in the barrier structure^ 9 -, which contribute collectively to the electric current 
and shot noise. Thus, the current is lowered at the peak and simultaneously lifted in both 
sides around the peak in the current-bias spectra. When the SOC is present, spin "f 
electrons exclude spin "J," electrons as well as spin ones. Therefore, the large noise 
suppressions are a consequence of the repulsion between current pulses of different spin 
states in addition to the consequence of the Pauli blockade and Coulomb repulsion. 

III. QUANTUM PUMPING BEYOND LINEAR RESPONSE 
A. Introduction to quantum pumping 

Generally speaking, the transport of matter from low potential to high potential excited 
by absorbing energy from the environment can be described as a pump process. The driving 
mechanics of classic pumps is straightforward and well understood 3 ^. The concept of a 
quantum pump is initiated several decades ago 3 ^ with its mechanism involving coherent 
tunneling and quantum interference. Research on quantum pumping has attracted heated 
interest since its experimental realization in an open quantum dot^~— . 

The mechanisms of an adiabatic quantum pump can be demonstrated in a mesoscopic 
system modulated by two oscillating barriers (see Fig. 9). To prominently picture the charge 
flow driven process within a cyclic period, the two potential barriers are modulated with a 
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phase difference of vr/2 in the manner of XJ\ = U + Ui u sint and U 2 = U + U 2u sin(t + vr/2). 
Our discussion is within the framework of the single electron approximation and coherent 
tunneling. The Pauli principle is taken into account throughout the pumping process. The 
Fermi energy of the two reservoirs and the inner single-particle state energy are equalized 
to eliminate the external bias and secure energy-conserved tunneling [The kinetic properties 
(charge current, heat current, etc.) depend on the values of the scattering matrix within the 
energy interval of the order of max(fcgT, fruS) near the Fermi energy. In the low-frequency 
{oj — > 0) and low-temperature (T — > 0) limit we assume the scattering matrix to be energy 
independent]. As shown in Fig. 9, the transmission strengths between one of the reservoirs 
and the inner single-particle state are denoted by ti-t4. When t G [0,7r/2], sint changes 
from to 1 and sin(t + vr/2) changes from 1 to 0. Considering the time-averaged effect, the 
chance of U\ > U 2 and U\ < U 2 is equal. Therefore, the probability of t 1 and t 3 balance 
out. The tunneling quantified by t 2 and £4 do not occur since the inner particle state is 
not occupied. When t G [vr/2,vr], sint changes from 1 to and sin(£ + vr/2) changes from 
to -1. U\ > U 2 invariably holds in this time regime. The probability of £3 prevails and a 
net particle flow is driven from the right reservoir to the middle state. When t G [vr,3vr/2], 
sint changes from to -1 and sin(t + vr/2) changes from -1 to 0. The probability of t 2 
and £4 balance out and the tunneling quantified by t\ and £3 are excluded from the Pauli 
principle. No net time-averaged tunneling occurs. When t G [3vr/2,27r], sint changes from 
-1 to and sin(t + vr/2) changes from to 1. U\ maintains a lower height than U 2 , which 
drives the particle in the inner state to the left reservoir. Through one whole pumping 
cycle, electrons are pumped from the right reservoir to the left by absorbing energy from 
the two oscillating sources. The tunneling is governed by quantum coherence. In each 
period, the pumping process repeats and the particles are driven continuously in the same 
direction as time accumulates. Direction-reversed pumped current can be obtained with 
reversed phase difference of the two oscillating gates. The direction of the pumped current 
is from the phase-leading gate to the phase-lagged one without exception when we assume 
that higher barriers admit smaller transmission probability. It can find resemblance in its 
classical turnstile counterpart^ in which the fore-opened gate admits transmission ahead of 
the later-opened one driving currents in corresponding manner. 

The current and noise properties in various quantum pump structures and devices 
were investigated such as the magnetic-barrier-modulated two dimensional electron gas^, 
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mesoscopic one-dimensional wire^ 1 ^, quantum-dot structures^s^ 1 ^ 1 ^, mesoscopic rings 
with Aharonov-Casher and Aharonov-Bohm effect^, magnetic tunnel junctions^, chains 
of tunnel-coupled metallic islands^, the nanoscale helical wire^,the Tomonaga-Luttinger 
liquid 59 , and garphene-based devices^ 5 - 1 ^. Theory also predicts that charge can be pumped 
by oscillating one parameter in particular quantum configurations 58 . A recent experiment^ 
based on two parallel quantized charge pumps offers a way forward to the potential appli- 
cation of quantum pumping in quantum information processing, the generation of single 
photons in pairs and bunches, neural networking, and the development of a quantum stan- 
dard for electrical current. Correspondingly, theoretical techniques have been put forward 
for the treatment of the quantum pumps^^ 5 ^ 1 ^ 1 ^ 1 ^. One of the most prominent is the 
scattering approach proposed by Brouwer who presented a formula that relates the pumped 
current to the parametric derivatives of the scattering matrix of the system. Driven by 
adiabatic and weak modulation (the ac driving amplitude is small compared to the static 
potential), the pumped current was found to vary in a sinusoidal manner as a function of 
the phase difference between the two oscillating potentials. It increases linearly with the 
frequency in line with experimental finding. The Floquet scattering theory is developed^ 
for quantum-mechanical pumping in mesoscopic conductors. It can be used to investigate 
quantum pumping behavior at arbitrary pumping amplitude and frequency. 

As an example to demonstrate the Floquet scattering theory, we focus on the experimen- 
tally observed deviation from the weak-pumping theory with only the first-order parametric 
derivative of the scattering matrix considered. By expanding the scattering matrix to higher 
orders of the time and modulation amplitude, experimental observation can be interpreted 
by multi-energy-quantum-related processes. 

B. Theoretical formulation 

We use the scattering matrix approach to describe the response of a mesoscopic phase- 
coherent sample to two slowly oscillating (with a frequency u) external real parameters Xj(t) 
(gate potential, magnetic flux, etc.), 

Xj (t) = X 0>j + X^-e^-w) + X^e-^-w), j = i, 2 . (31) 
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Xqj and X u j measure the static magnitude and ac driving amplitude of the two parameters, 
respectively. The phase difference between the two drivers is defined as = (p 1 — ip 2 . The 
mesoscopic conductor is connected to two reservoirs at zero bias. The scattering matrix s 
being a function of parameters Xj(t) depends on time. 

A time-dependent scattering matrix can be introduced as follows. 

/oo 
S a p(t,t')ap(t')dt', t>f. (32) 
-oo 

Its Wigner transform reads 

/oo 
e iE ^S aP (t : t')dt'. (33) 
-oo 

We assume the scattering time t — t' is small. Up to corrections of order ftu/^ (7 measures 
the escape rate), the matrix S a p{E,t) is equal to the "instantaneous" scattering matrix 
S X (E), which is obtained by "freezing" all parameters Xj to their values at time t. Below, 
we use the instant scattering matrix s (t) in place of S a p (E, t) to describe the physics for 
simplicity. The kinetic properties (charge current, heat current, etc.) depend on the values 
of the scattering matrix within the energy interval of the order of max(/csT, t\w) near the 
Fermi energy. In the low-frequency (uu — > 0) and low-temperature (T — > 0) limit we assume 
the scattering matrix to be energy independent. To investigate the deviation from the small 
amplitude X u j limit, we expand the scattering matrix s(t) into Taylor series of Xj(t) to 
second order at X j with the terms linear and quadratic of X^j present in the expansion, 

s (t) « s (X 0J ) + L u e ,wt + s+ve-™ + s 2 + s^e 2 ^ + s +2u] e- 2 ^\ (34) 

with 

' s ±UJ = £ X^^ds/dXj, 
3=1,2 

< S2 = E X^dH/dXl (35) 
3=1,2 

s± 2 . = \ E xl^dH/dx 2 . 

V 3=1,2 

It can be seen from the equations that higher orders of the Fourier spectra enter into the 
scattering matrix. As a result, both the nearest and next nearest sidebands are taken into 
account, which implies that a scattered electron can absorb or emit an energy quantum of 
hu or 2hu before it leaves the scattering region. In principle, third or higher orders in the 
Taylor series can be obtained accordingly. However, the higher-order parametric derivatives 
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of the scatter matrix diminish dramatically and approximate zero. Numerical calculation 
demonstrates that even in relatively large amplitude modulation, their contribution is neg- 
ligible. 

The pumped current depends on the values of the scattering matrix within the energy 
interval of the order of max (ksT, 2hu) near the Fermi energy. In the low-temperature limit 
(T — > 0), an energy interval of 2ftw is opened during the scattering process. 

The mesoscopic scatterer is coupled to two reservoirs with the same temperatures T 
and electrochemical potentials /i. Electrons with the energy E entering the scatterer are 
described by the Fermi distribution function fo(E), which approximates a step function at 
a low temperature. Due to the interaction with an oscillating scatterer, an electron can 
absorb or emit energy quanta that changes the distribution function. A single transverse 
channel in one of the leads is considered. Applying the hypothesis of an instant scattering, 
the scattering matrix connecting the incoming and outgoing states can be written as 

ba(t)=^2s aP (t)a p (t). (36) 

P 

Here s a p is an element of the scattering matrix s; the time-dependent operator is a a (t) = 
J dEa a (E) e - lEt / h ^ and the energy-dependent operator a a (E) annihilates particles with to- 
tal energy E incident from the a lead into the scatter and obey the following anticommutation 
relations 

[at (E) , a p (E')] = 5 a ?5 (E - E') . (37) 

Note that above expressions correspond to single- (transverse) channel leads and spinless 
electrons. For the case of many-channel leads each lead index (a, /3, etc.) includes a 
transverse channel index and any repeating lead index implies implicitly a summation over 
all the transverse channels in the lead. Similarly an electron spin can be taken into account. 
Using Eqs. (34) and (36) and after a Fourier transformation we obtain 

b a (E) = [h,<xpap (E) + S2, a pap (E) + S-^fi/j (E + hu) 

(38) 

+s +0J , a pap (E-hu}) + s-2u,apap (E + 2hu) + s +2Ut apap (E - 2fiw)] . 

The distribution function for electrons leaving the scatterer through the lead a is fa^ (E) = 
b\ x {E)b a (E)\ where (••■) means quantum-mechanical averaging. Substituting Eq. (38) 
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we find 

jr(out) ^ _ j2 [|s 0jQ( g + S 2 ,al3\ 2 fo (E) + |s_^ Q/3 | 2 f (E + ftw) 

(39) 

I^+oj, a/3 

\ 2 f (E + 2hw) + \s | 2 /o(B-2M]- 
The distribution function for outgoing carriers is a nonequilibrium distribution function 
generated by the nonstationary scatterer. The Fourier amplitudes of the scattering matrix 
|s-o;,a/3| 2 (|s+oj,a^| 2 ) is the probability for an electron entering the scatterer through the lead 
(3 and leaving the scatterer through the lead a to emit (to absorb) an energy quantum 
hu and \s>-2uj,ap\ 2 (|s+2w,a,s| 2 ) is that of the energy quantum 2ftoj process. \§o, a p + S2,a,s| 2 is 
the probability for the same scattering without the change of an energy with the second- 
order term <§2,a/3 much smaller than the zero-order term So,a/3 in weak-modulation limit 
(X u j <C Xqj) and can be omitted therein. 

Using the distribution functions fo{E) for incoming electrons and f° ut (E) for outgoing 
electrons, the pumped current measured at lead a reads 

^ = ^1°° (ft ( E ) h * (£)) " < a « ^ W) dK ( 4 °) 
Substituting Eqs. (34) and (30) we get 

T ew_ ST V V ® s a0 ® S a0 Q • • / \ 

1 P — 2tt ^uji-^-Ujh dXj, 8Xj„ Yh) 

0,jl ' h * av ( 41 ) 

+f £ XZ^AJ&H&i ^ [2 (^ - <p h )]. 

P,hd2 31 32 

Quantum pumping properties beyond the theory based on first-order parametric derivative 
of the scattering matrix are demonstrated in Eq. (41). By taking higher orders of the Fourier 
spectrum of the scattering matrix into consideration, double ftw energy quantum (or a 2hu 
energy quantum) emission (absorption) processes coact with single Huo quantum processes. 
In the weak-modulation limit, the second term in the right-hand side of Eq. (41) is small, 
which implies that double hw quantum processes are weak and therefore not observable. As 
the ac driving amplitude is enlarged, this term increases markedly and contribution from 
double hu quantum processes takes effect. As a result, the dependence of the pumped 
current on the phase difference between two driving oscillations deviates from sinusoidal 
and changes from sin to sin 20, which is observed in experiment^. Moreover, the relation 
between the pumped current and the ac driving amplitude X u j is reshaped. It is also seen 
that the linear dependence of the pumped current on the oscillation frequency holds for 
multi-quanta-related processes. 
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C. Numerical results and interpretations 

Here, numerical results of the pumped current in a two-oscillating-potential-barrier mod- 
ulated nanowire are presented and comparison with experiment is given. We consider a 
nanowire modulated by two gate potential barriers with equal width L = 20 A separated by 
a 2L = 40 A width well (see Fig. 10). The electrochemical potential of the two reservoirs 
H is set to be 60 meV according to the resonant level within the double-barrier structure. 
The two oscillating parameters in Eq. (31) correspond to the two ac driven potential gates 
-Xi,2 if) U± t 2 (t) with all the other notations correspond accordingly. We set the static 
magnitude of the two gate potentials Uo,t = t/0,2 — Uq — 100 meV and the ac driving 
amplitude of the modulations equal U U) i = U Wj 2 = U w . 

In Fig. 11, the dependence of the pumped current on the phase difference between the 
two ac oscillations is presented. In weak- modulation regime (namely U u is small), sinusoidal 
behavior dominates. Here, three relatively large is selected to reveal the deviation from 
the sinusoidal dependence. (The magnitude of the pumped current mounts up in power- 
law relation as a function of as shown in Fig. 12. The sinusoidal curve for small U u 
would be flat and invisible in the same coordinate range.) It can be seen from the figure 
that the J p -0 relation varies from sinusoidal (sin0) to double-sinusoidal (sin 20) as the ac 
oscillation amplitude is increased. The interpretation follows from Eq. (41). The single hw 
quantum emission (absortion) processes feature a sinusoidal behavior while the 2hw quantum 
emission (absortion) processes feature a double-sinusoidal behavior when the Fourier index is 
doubled. As U u is increased, double hw quantum processes gradually parallel and outweigh 
the single hw quantum ones. It is also demonstrated that when the single hw quantum 
processes have the effect of sin0 dependence, the double hw quantum processes induce a 
— sin 20 contribution with a sign flip, which can be understood from the sign change of the 
derivative of the scattering matrix. The effect of three- and higher hw quantum processes is 
small even for large U u comparable to Uq. The experimental observations^ as a deviation 
from the weak-modulation limit are revealed by our theory. 

Experiment^ also discovered that for weak pumping the dependence of the pumped 
current on the pumping strength obeys a power of 2 relation, as expected from the simple 
loop-area argument^; for strong pumping, power of 1 and 1/2 relation is observed. We 
presented in Fig. 12 the numerical results based on our theory of the I p -U u relation at a 
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fixed 0. To demonstrate its power-law dependence, natural logarithm of the variables is 
applied. From Eq. (41), it can be seen that for large ac driving amplitude U w , contribution 
of double fajj quantum processes (formulated in the second term on the right hand side 
of the equation) causes the I p -U u relation to deviate from its weak-modulation limit, the 
latter of which is I p oc U%. For different phase difference between the two ac drivers, the 
deviation is different. At = ir the pumped current is invariably zero regardless of the order 
of approximation determined by time-reversal symmetry. At = n/2, sin 20 is exact zero, 
and no difference is incurred by introducing higher order effect. If we shift the value of to 
0.497T, the abating effect of the double Tujj quantum processes has the order of U* with the 
small second-order parametric derivative of the scattering matrix smoothing that effect a 
bit. Consequently, a power of 2 — > 1 — > 1/2 relation is obtained and visualized by the curve 
fit, which is analogous to experimental findings. For different values of 0, sharper abating 
and augmental effect occurs with analogous mechanisms. It is possible that the experiment^ 
was done at the phase difference close to tt/2 while trying to approach maximal pumped 
current in the adiabatic and weak-pumping limit. 



IV. SUMMARY AND FUTURE DIRECTIONS 



The scattering matrix method is initiated by Landauer and Biittiker to investigate the 
conductance of multi-terminal and multi-channel mesoscopic conductors. The spin degree of 
freedom can be included in the formalism by enlargement of the dimension of the scattering 
matrix. The current-current correlation and spin-spin correlation, such as the shot noise, can 
be calculated from the cross products of the scattering matrix. Along this direction, higher- 
order correlation function can also be considered. Dynamic transport processes including the 
quantum pumping behavior can be dealt with by the time-dependent scattering approach. 

The development of the scattering theory enables its potential applications in currently 
open issues. The interaction can be included to the scattering matrix by the renormalization 
factors. The shot noise properties in various conductors with active spin degree of freedom 
can be considered. The time-dependent scattering theory provides a way to deal with dy- 
namic quantum issues. Some particular problems include non-harmonically driven quantum 
pumping, spin pumping in racetrack memory applications, and multiferroic transport dy- 
namics, etc.. 
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FIG. 1: Schematic demonstration of a single-barrier tunneling problem. The quantum states 
between two reservoirs at zero bias is indicated. 
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FIG. 2: Schematic demonstration of a quantum wire modulated by two potential barriers. The 
quantum state between two reservoirs at zero bias is indicated. 
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FIG. 3: Time-integrated current density as a function of the phase difference between the two 
modulations for different modulation frequencies. The Fermi level of the two reservoirs Ep = 60 
meV counting from the conduction band edge of the electron gas structure. 
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FIG. 4: Time variation of the current density within an oscillating cycle. Different phase difference 
between the two modulations is considered. The modulation frequency is set to be w = 10 MHz 
and the Fermi level Ep to be 60 meV. 
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FIG. 5: Schematics of the double-barrier resonant diode. The resonant level is sketched between 
the two barriers. The upper panel demonstrates the resonant level in conventional diode without 
the spin-orbit coupling (SOC). In the lower panel, the SOC behaves like a pseudomagnetic field 
and induces a split of different spin components of the resonant level in the barrier structure, which 
contribute collectively to the electric current and shot noise. 



31 



L 




FIG. 6: Schematics of the scattering approach. Spin components of the incoming and outgoing 
states are indicated. 
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FIG. 7: Current / (a) and shot noise S (b) as functions of the applied bias EV of electrons 
traversing a symmetric DBRD structure with different Dresselhaus constants 7. 
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FIG. 8: Fano factor as a function of the applied bias EV (a) and S /2eI norm VS I '/Inarm (b) of 
electrons traversing a symmetric DBRD structure with different Dresselhaus constants 7. Two 
straight lines in (b) show the full shot noise value (2eJ) and half of its value for comparison. 
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FIG. 9: The tunneling scenario of an adiabatic quantum pump. The two shadowed blocks represent 
the left and right electron reservoirs respectively. The two barriers oscillate adiabatically in time. 
The middle bar indicates the single-particle state between the two barriers. The Fermi levels of 
the two reservoirs are the same and are leveled to the single-particle state within the conductor. 
t±-t4 indicate the transmission amplitudes between one of the two reservoirs and the middle single- 
particle state. 



35 




FIG. 10: Schematics of the quantum pump: a nanowire modulated by two ac driven potential 
barriers. 
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FIG. 11: Pumped current as a function of the phase difference between the two modulations for 
different ac driving amplitudes. 
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FIG. 12: Pumped current as a function of the ac driving amplitude U u along with fits to I p oc 
U^j (red solid circle) below 35 meV, I p oc U w (green upward triangle) below 41 meV, and I p oc 

1/2 

Uj above 41 meV (blue downward triangle). To demonstrate its power- law dependence, natural 
logarithm of the variables is applied. The phase difference between the two ac driver <f> = 0.497T. 
Inset is the zoom-in of the circled region. 
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